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Dyson orbitals play an important role in understanding quasi-particle effects in the correlated 
ground state of a many-particle system and are relevant for describing the Compton scattering cross 
section beyond the frameworks of the impulse approximation (1A) and the independent particle 
model (IPM) . Here we discuss corrections to the Kohn-Sham energies due to quasi-particle effects in 
terms of Dyson orbitals and obtain a relatively simple local form of the exchange-correlation energy. 
Illustrative examples are presented to show the usefulness of our scheme. 
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I. INTRODUCTION 

Dyson orbitals are a set of one-particle orbitals that are 
associated with many-electron wavefunctions. These or- 
bitals connect the exact ground-state of the JV-electron 
system with excited states containing TV — 1 or N + 1 
electrons. The importance of Dyson orbitals in under- 
standing Compton scattering spectra has been empha- 
sized recently by Kaplan et al. |l| , who present a general 
formalism for the Compton cross section, which goes be- 
yond the standard treatment involving the frameworks 
of the impulse approximation (IA) and the independent 
particle model (IPM). The breakdown of the IA in de- 
scribing core Compton profiles is well-known |2|. More 
recently, high resolution valence Compton profiles (CPs) 
of Li at relatively low photon energy of 8 — 9 kcV have 
been found to show asymmetries in shape and smearing 
of the Fermi surface (FS) features where deviations from 
the I A have been implicated 0, 0] . 

Dyson orbitals also give insight into quasi-particle ef- 
fects in the correlated ground state of the many body 
system. In this article, we focus on understanding ener- 
gies of one-particle excitations of the ground state, which 
play an important role not only in the formalism of the 
Compton scattering cross-section, but also in the band 
structure problem more generally. To this end, a Green's 
function approach is used to first obtain an expression 
for the exchange-correlation energy in terms of the self- 
energy operator involved in the description of the Dyson 
orbitals. A local ansatz for the self-energy is then invoked 
to obtain a relatively simple expression for the excita- 
tion energies. We illustrate our scheme by considering 
the example of first ionization energies of low Z atoms 
from Z = 1 (H) to Z — 6 (C) and find good agreement 
with the corresponding experimental results. As another 
example, the measured bandgap in diamond is also re- 
produced reasonably by our computations. 

An outline of this article is as follows. The introduc- 
tory remarks are followed in Section II by a brief overview 
of the general formalism of Ref . for the Compton scat- 
tering cross section. The importance of properly includ- 
ing excitation energies in the computation for describing 



the asymmetry of the CP around q = is stressed. Sec- 
tion III presents the Green's function formulation and 
discusses our scheme for computing excitation energies. 
Section IV gives a few illustrative applications of the the- 
ory, followed in Section V by a few concluding remarks. 

II. DYSON ORBITALS AND COMPTON 
SCATTERING 

The Dyson spin-orbital g n can be defined in terms of 
the many-body ground-state wavefunction an d the 
wavefunction ty n of the singly ionized excited system 
characterized by the quantum number n ]±&, LI LD, M, ■ 

g n (*N) = y/N J *„(xi....xjv-i)**o(xi....xat) dx!....dx N 

(1) 

where the index x denotes both spatial and spin coor- 
dinates and the integration over dni implicitly includes 
a summation over the spin coordinates. The Dyson or- 
bitals thus give generalized overlap amplitudes between 
the ground state and the singly ionized states of the many 
body system. Note that, in general, Dyson orbitals do 
not form an orthonormal set. The Dyson spin-orbital 
with the spin projection a may be written in terms of 
the spin function cr(C) as 

<7„(x) =ff n (r,<7(0) =<?„(r)a(C) . (2) 

The excitation energy E^ associated with the n th Dyson 
orbital is given by 

= E n (N - 1) - E Q (N) (3) 

where Eq(N) is the Af-particle ground state energy and 
E n (N — 1) is the energy of the (N — 1) particle ionized 
system when it is in its n th quantum state. Eqs. 1-3 
above define what may be thought of as occupied Dyson 
orbitals. For completeness, one can also introduce a par- 
allel set of "unoccupied" Dyson orbitals, which connect 
the ground state to various (N + 1) particle states con- 
taining an added electron. 

The importance of Dyson orbitals in understanding 
the nature of the Compton scattering spectra beyond the 



2 



framework of the IA and the IPM has been emphasized 
by Kaplan et al. Q. In particular, the triple-differential 
scattering cross-section for the (7,67) process, which is 
the elementary process underlying Compton scattering, 
is given by 



d 3 G 



duj2(lVL2dVL 
^2 lffn(q)| 2 ^l —V2 — E, 
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(4) 



where gn(q) is the Fourier transform of g n (r), q is the 
momentum transfered to the final system, u>i and 0J2 are 
respectively the energies of the photon before and after 
the collision, and the summation extends over the oc- 
cupied Dyson orbitals. Eq. (@J assumes a large energy 
transfer (u>i — u>2» E^), so that the outgoing electron 
possesses sufficiently high energy that it can be approx- 
imated by a free electron plane wave form. The bind- 
ing energy E b in the (5-function in Eq. Q is usually 
neglected in obtaining the standard IA-based expression 
for the cross section. The resulting form of the Compton 
profile after coordinates of the recoil electron are inte- 
grated over can be shown to be symmetric around q = 0. 
Therefore, an accurate computation of the binding ener- 
gies E^ is important for describing the asymmetry of 
the Compton profile related to deviations from the IA. 



III. GREEN'S FUNCTION METHOD 

In order to gain a handle on the nature of the excitation 
energies E^\ it proves useful to approach the problem 
through the Green's function method. The orbitals <7 n (r) 
of Eq. 0) satisfy the Dyson equation [3, 0, E| 



(^+V ext (r)+V H (r))g n (r) 
+ J dV E xc {r,r',E^)g n (r') = E^g n {v), (5) 

where V ex t (r) is the external potential, Vh is the Hartree 
potential and Yi xc is the self-energy. 

The Green's function can be expressed as 



G(r,r»=£ 



co - Ei n) + zSsign(E { b n) - p) ' 



(6) 



where p is the chemical potential and 6 is an infinitesimal 
positive number. The density matrix p(r, r') and the 
electron density distribution n(r) = p(r, r) are obtained 
by integrating the spectral function A(r,r',uj) over the 
occupied electronic energies 



p{r,r')= f duA(r,r',oj) 

J —OO 



(7) 



where the spectral function A(r, r',w) is given by 

A(r,r» =^ 9n (r)s;(r')- - 

ir[{uj-E y b 'Y + S 2 } 



(8) 



The standard calculation of the self-energy H xc and 
the Green's function G proceeds via the many-body per- 
turbation theory (MBPT) [T^]. The first order in the 
MBPT leads to the so-called GW approximation 
(G stands for the Green's function and W denotes the 
screened Coulomb interaction). The GW equation for 
the self-energy is 



£ xc = iG W , 



(9) 



where Go is the th order Green's function. The cal- 
culation of W usually requires heavy computations of 
the dielectric function. Here, we will take a shortcut by 
describing the screening in terms of the pair-correlation 
function. 

Since many-body effects beyond the Hartree approxi- 
mation are contained in the exchange-correlation energy 
E xc , we can write 



1 

E xc = - / dujd 3 rd 3 r'T, xc (r, r', oj)A(r, r', w) 

" J -co 



(10) 



where the integral is over the occupied states. In the 
Hartree-Fock limit, the self energy operator is instante- 
nous, i.e., E xc (r,r',oj) — S x (r, r'). By integrating over 
to, ones finds the well known result for the exchange en- 
ergy E x 



E X = \J d 3 rd 3 r'Z x (r,r')p(v,r') 



(11) 



We now consider the standard expression of the 
exchange-correlation energy E xc in terms of the 
exchange-correlation hole |15| 



Err = - I d 3 rd 3 r 



1 dX n(r)n(r')C x (r,r>) 



(12) 



where A is the coupling constant from the Hellmann- 
Feynman theorem, and C\ is a pair-correlation function 
describing the exchange-correlation hole. The two ex- 
pressions of Eqs. (|10[1 and l|12Jl for E xc can be linked by 
assuming a local and instantaneous ansatz for the self- 
energy operator 



£ xc (r,r» = 2e xc (r)<5(r - r') 



Here 



eUr)= lS drl L 



d\ 



n(r)n(r')C A (r, r') 



(13) 



(14) 



is the exchange-correlation energy per particle. 

Ea. i|13fl has an intuitive interpretation. When an elec- 
tron moves, its exchange-correlation hole moves with it 
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and modifies its effective mass and excitation energy. In- 
sight into the local approximation for Y, xc is provided by 
considering the homogeneous electron gas. In this case, 
one obtains 



E xc (r,r' ) w) = 2e xc £(r-r'), 



where e xc is a constant due to the translation invariancc 
of the system. A simple expression for e xc in the metallic 
density range is given by [17J 



0.916 0.127 



Ry, 



(16) 



where the mean interelectronic spacing r s is obtained via 
the electron density n by n(47rri?/3) = 1. Thus, as a 
consequence of the short range of the self-energy, the 
electron-electron interaction leads only to a uniform shift 
of the levels in relation to the non-interacting gas. In the 
Hartree-Fock case, the density of states unphysically goes 
to zero at the Fermi level as a result of the long range of 
the self-energy The situation in a real metal is pre- 
sumably closer to that of the homogeneous electron gas 
as the screening tends to make the correct T, xc local. In- 
deed, band calculations based on the Density Functional 
Theory (DFT) employing local potentials [l5j have been 
rather successful in reproducing the experimental FS's in 
wide classes of materials [lij ]. 

Bearing these considerations in mind, the correction to 
the excitation energy of Eq. © may be obtained in 
the first perturbational order as 

E { b n) = e n + J d 3 r (2e xc (r) - v xc (r))\g n (r)\ 2 , (17) 

where e n denotes the Kohn-Sham eigenvalue and 
Vxci?) = SE XC /Sn(r) is the exchange-correlation poten- 
tial in the Kohn-Sham equations. (/ n (r) on the right hand 
side of Eq. (|17|l can be reasonably replaced by the Kohn- 
Sham orbitals [IJj . The exchange-correlation energy per 
particle e xc is computed straightforwardly within the lo- 
cal density approximation (LDA) 0. Notably, correc- 
tion to E^ in Eq. (|17ll will in general be state-dependent 
and therefore this spectrum will be a useful starting 
point for implementing the scheme proposed by Barbi- 
ellini and Bansil 20] for treating the momentum density 
and Compton profile of the correlated ground state of 
the anisotropic electron gas beyond the Lam-Platzman 
correction. 



IV. EXCITATION ENERGY CALCULATIONS 

Fig. ^ provides an illustrative example of the useful- 
ness of Eq. (|17|l . The first ionization energy of atoms 
from H to C (Z = 1 — 6) is considered using rela- 
tivists DFT atomic wavefunctions [2l|. The exchange- 
correlation energy e xc (r) and the potential v xc (r) have 



(is) r 



>. 

c 

uj ! 

c 
o 

(0 
N 

"E 
o 

- 0.5 












Experiment 








O LDA 










+ QP correction 


- 


o 






- 













+ 








* 








V 




o 






o 


O 






o 




O 



1 2 3 4 5 

Atomic number (Z) 



FIG. 1: Comparison of the theoretical and experimental 
atomic first ionization energies for H {Z = 1) to C (Z = 6). 
Diamonds are the experimental values |23|. open circles give 
LDA values, while crosses give values corrected by using Eq. 
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been calculated within the LDA parametrized by Hcdin 
and Lundqvist 22] . The LDA eigenvalues (open circles) 
are seen to be substantially lower than the experimental 
values (diamonds). The quasiparticle correction of Eq. 
(|17fl brings the theoretical values (crosses) in substantial 
agreement with the experimental values 23]. Note that 
in the case of H, although the quasiparticle correction 
yields a large improvement, the self-interaction error is 
not fully removed and about 20 % discrepancy between 
theory and experiment remains. 

We have also applied our scheme to investigate the 
band gap in diamond. For this purpose, the selfconsistent 
electronic band structure of diamond |24j was obtained 
using the local density exchange-correlation functional of 
Hedin and Lundqvist As expected, the LDA band 

gap of 3.99 eV so computed is too small. The inclusion 
of the correction of Eq. I|17|l yields a gap of 5.3 eV, 
which is in far better agreement with the experimental 
value of 5.48 eV j^. More generally, we expect our 
correction to reproduce semiconductor bandgaps with an 
accuracy comparable to that of the com put ationally more 
involved GW approximation [2^E7lE3 l29 l l3C | . It should 
also be noted that in analyzing the generalized density- 
functional theory (GDFT), several authors have pointed 
out that the energy of electronic excitations across the 
gap of insulators and semiconductors can be expressed as 
the sum of the so-called Kohn-Sham gap and a correction 
that is usually of the same order of magnitude [3lll32ll33] ]. 

We note that quasi-particle effects directly influence 
the shape of the FS of a metal. For instance, the FS of 
Cu can be measured very precisely so that direct com- 
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parisons with theoretical predictions are easy to follow, 
de Hass-van Alphen (dHvA) measurements |3J] find the 
large belly of the FS of Cu to be more spherical compared 
to the LDA computations. In V, dHvA [3{| and positron 
annihilation pjjj experiments indicate that the FS pock- 
ets at the N symmetry point are smaller than those from 
LDA calculations. In 7-Ce, an improved agreement with 
positron annihilation results is obtained with the LDA 
FS if the /-band is moved up by 40 — 50 mRy 

H3. In- 

terestingly, aforementioned discrepancies concerning the 
FSs can be substantially corrected via a simplified com- 
putation of the self-energy via the exchange correlation 
hole !H. 



V. SUMMARY AND CONCLUSIONS 

We discuss aspects of Dyson orbitals for gaining insight 
into quasi-particle effects in the correlated ground state 
of a many-particle system. The importance of Dyson or- 
bitals, which connect the many-body ground state with 
its singly ionized excited states, has been emphasized pre- 
viously for describing Compton scattering profiles beyond 
the limitations of the IA and the IPM and we start with 
a brief review of this earlier study 0] . We focus on de- 
lineating corrections to the Kohn-Sham energies due to 
quasi-particle effects and show how one can gain a han- 
dle on obtaining substantially improved excitation ener- 
gies in molecules and solids. For this purpose, a Green's 
function approach is used to suggest a relatively simple 
local form of the exchange-correlation energy. Illustrative 
applications of our scheme indicate that our self-energy 
expression should provide a reasonable description of ex- 
citation energies in metals and semiconductors with an 
accuracy roughly comparable to that of the GW method. 

This work is supported by the US Department of En- 
ergy contract DE-AC03-76SF00098 and benefited from 
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Northeastern University's Advanced Scientific Computa- 
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